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Abstract 

This note addresses the homogenization error for linear elliptic equations in divergence- 
form with random stationary coefficients. The homogenization error is measured by com¬ 
paring the quenched Green’s function to the Green’s function belonging to the homogenized 
coefficients, more precisely, by the (relative) spatial decay rate of the difference of their 
second mixed derivatives. The contribution of this note is purely deterministic: It uses 
the expanded notion of corrector, namely the couple of scalar and vector potentials r), 
and shows that the rate of sublinear growth of (0, a) at the points of interest translates 
one-to-one into the decay rate. 


1 A brief overview of stochastic homogenization, and a common 
vision for quenched and thermal noise 

Heterogeneous materials typically have the behavior of a homogeneous material on scales that 
are large with respect to the characteristic length scale of the heterogeneity, provided that the 
“type of heterogeneity” is the same all over the medium. Here, we think of material properties 
like conductivity or elasticity that are typically described by an elliptic differential operator in 
divergence form —V • aV, where the heterogeneity resides in the uniformly elliptic coefficient 
field a. Homogenization then refers to the fact that on large scales, the solution operator 
corresponding to —V • aV behaves like to the solution operator to —V • a^V with a constant 
coefficient ah- 

This type of homogenization is well-understood in case of a periodic medium, where, for instance, 
ah can be inferred from a cell problem and the homogenization error can be expanded to any 
order in the ratio of the period to the macroscopic scale of interest. However, the case where only 
the statistical specification of the coefficient field is known might be more relevant in practice. 
In this case, one thinks of an ensemble (■) of coefficient fields and one speaks of stochastic 
homogenization. For homogenization to occur, the statistics have to be translation invariant, 
meaning that the joint distribution of the shifted coefficient field is identical to that of the 
original coefficient field, a property called stationarity in the parlance of stochastic processes. 
For homogenization to be effective, the statistics of the coefficient field have to decorrelate (more 
precisely, become independent) over large distances. Hence stochastic homogenization relies on 
the separation of scales between the correlation length and the macroscopic scale of interest. 

From a qualitative point of view, stochastic homogenization has been rigorously understood 
since the seminal works of Kozlov m and of Papanicolaou & Varadhan [21]: Stochastic ho- 
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mogenization takes place when the ensemble is stationary and ergodic, the latter being a purely 
qualitative way of imposing decorrelation over large distances. 

Stochastic analysis has a finer, still qualitative view, on stochastic homogenization, where 
—V • aV is seen as the generator of a (reversible) diffusion process. In a discrete medium 
this leads to the picture of a “random walk in random environment”, which amounts to the 
superposition of thermal noise (the random walk) and quenched noise (the random jump rates). 
Here, the relevant qualitative question is that of a “quenched invariance principle”: On large 
scales and for large times (on a parabolic scale), the random walker behaves like a Brownian 
motion (with covariance given by ah) for almost every realization of the random environment. 
Surprisingly, first full results in that direction came quite a bit later [22]. Stochastic analysis 
is still mostly interested in qualitative results, but pushing the frontier in terms of models, for 
instance towards degenerate diffusions like random walks on (supercritical) percolation clusters 
or towards diffusions with random drifts leading to non-reversible random walks. 

Numerical analysis has another, naturally more quantitative view on stochastic homogenization. 
As opposed to periodic homogenization, there is no cell problem to extract the homogenized 
coefficient ah- Hence one has to resort to an artificial “representative volume”: On such a cube 
(let us adopt three-dimensional language), one samples a realization of the medium according to 
the given statistical specifications and then solves three boundary value problems for —V • aS7u , 
corresponding to describing different slopes on u. In the case of a conducting medium, where 
the coefficient field a corresponds to the heterogeneous conductance, u corresponds to the elec¬ 
trical potential, so that the boundary conditions impose an average potential gradient, that is, 
an electrical field in one of the three coordinate directions. One then monitors the average elec¬ 
trical current aVu; this (linear) relation between average electrical field and average electrical 
current yields an approximation to the homogenized coefficient ah as a linear map. Clearly, this 
approximate homogenized coefficient ah l will be the closer to the true homogenized coefficient 
ah, the larger the linear size L of the cube is, where the relevant (small) nondimensional param¬ 
eter is the ratio of correlation length to L. Intuitively, there are two error sources: On the one 
hand, ah.L is a random quantity, since it depends on the given realization a of the coefficient 
field on the cube, so that there is a random error coming from the fluctuations of ah.L which can 
be reduced by repeated sampling. On the other hand, the very concept of the representative 
volume element perturbs the statistics, for instance in case of periodic boundary conditions, the 
concept introduces artificial long-range correlations, which is not affected by repeated sampling. 
Clearly, it is of interest to understand — which for mathematicians means to prove — how both 
errors scale in L. 

This natural and very practical question turns out to be difficult to analyze (rigorously). Shortly 
after the qualitative theory was introduced, Yurinskii [23] produced the first quantitative re¬ 
sult motivated by the above questions. He in fact used ingredients from stochastic analysis 
(the picture of diffusions in a random medium to understand sensitivities) and from regularity 
theory (Nash’s upper heat kernel bounds), but only obtained sup-optimal results in terms of 
exponents. After using (qualitative) stochastic homogenization to understand the large-scale 
correlation structure of a statistical mechanics model (fluctuating surfaces) in [T9], Naddaf Sz 
Spencer [20] realized that tools from statistical mechanics (spectral gap estimate) can be used 
to obtain quantitative results in stochastic homogenization (for discrete media). With Meyers 
perturbation of the Calderon-Zygmund estimate they introduced a second tool from elliptic 
regularity theory into the field, which allowed them to obtain optimal variance estimates in 
the case of small ellipticity ratio (small contrast media); this type of result was subsequently 
extended by Conlon and coworkers, see e. g. [6]. In [13, [2], Gloria and the last author used 
the same tool from statistical mechanics but yet another ingredient from elliptic regularity 
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theory (Caccioppoli’s estimate to obtain optimal spatially averaged bounds on the gradient of 
the quenched Green’s function) to obtain the first optimal error estimate on the representa¬ 
tive volume method also for large ellipticity ratios. In [IS], Marahrens and the last author 
used the concentration of measure property coming from the logarithmic Sobolev inequality to 
study the (random part of the) homogenization error itself, in form of optimal estimates on the 
variance of the quenched Green’s function. Using Green’s function estimates and two-scale ex¬ 
pansion, Gloria, Neukamm, and the last author m compared the heterogeneous and corrected 
homogenenous solution. While in [10] the error is measured in H 1 and averaged both over 
the domain and the ensemble, Gu and Mourrat m recently combined probabilistic techniques 
with Green’s function estimates to obtain a pointwise bound for solutions of both elliptic and 
parabolic equations. Since then, there has been a flurry of activities, which will be partially 
addressed in the next two sections, with the work of Armstrong & Smart [2] playing a central 
role. For instance, by now it is already understood that the error in the representative volume 
method is to leading order Gaussian m- We do not even mention the numerous activities in 
stochastic homogenization of non-divergence form equations, like fully nonlinear equations or 
Hamilton-Jacobi equations. 

In both qualitative and quantitative homogenization, both from the PDE and the stochastic 
analysis point of view, the notion of the corrector, a random function fa for every coordinate di¬ 
rection % = 1,..., d, is central. There is a very geometric and deterministic view of the corrector: 
Given a realization a of the coefficient field, x i —> x\ + fa(x) provides a-harmonic coordinates. 
This is also its main merit from the almost-sure stochastic analysis point of view: Seen in these 
coordinates, the diffusion turns into a Martingale. The corrector is also natural from the nu¬ 
merical analysis point of view: In the representative volume method, one actually solves for an 
approximate version of the corrector. Last not least, in the original (very functional analytic) 
PDE approach to stochastic homogenization the corrector is central: Using stationarity, one lifts 
the equation for the corrector to the probability space, solves it by Lax-Milgram, and expresses 
ah in terms of it. Like the work [8] on higher-order Liouville principles, this note demonstrates 
the usefulness of the vector potential a of the corrector, an object known in periodic homoge¬ 
nization, and recently introduced into random homogenization by Gloria, Neukamm, and the 
last author in [9]. 

We have a common vision for a regularity theory of random elliptic operators as considered 
in this note and for stochastic partial differential equations (SPDE). In other words, we’d like 
to capitalize more on the similarities between quenched and thermal noise. At first glance, 
these problems seem very different since in the first case, one is interested in the emergence 
of a generic large-scale regularity due to cancellations, while in the second case one wants to 
preserve a small-scale regularity despite the rough forcing. However, in both cases, the key is 
to understand how sensitively the solution <j) of an elliptic or parabolic equation (nonlinear to 
be of interest in case of SPDE with additive noise while in case of stochastic homogenization, 
the interesting effect is already present for a linear equations) depends on the data, be it the 
coefficients in case of stochastic homogenization or the right hand side in case of a random 
forcing, in which case this derivative can be associated to the Malliavin derivative. In this 
sensitivity analysis, one typically has to control the size of a functional derivative, that is, the 
functional derivative of some nonlinear functional of the solution (a norm, say) with respect to 
the data, which in case of the data being the coefficients is a highly nonlinear mapping even 
for a linear equation. Hence for a given realization of the data, one has to control the norm 
of a linear form on infinitesimal perturbations of the data. Since one is dealing with random 
stationary data, the appropriate measure of the size of the infinitesimal perturbations is best 
captured by an L 2 -type norm, the specific structure of which depends on the assumption on the 
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noise: An ordinary L 2 -norm in case of white-noise forcing or a more nonlocal (and thus weaker) 
norm in case of stochastic homogenization if one wants to cover also cases where the covariances 
of the coefficient field have a slow (that is, non-integrable) decay, as used in [5]. Even if the 
problem is a nonlinear one, the sensitivity estimates require a priori estimates for linear elliptic 
or parabolic equations, albeit with non-constant coefficients that a priori are just uniformly 
elliptic — it is here where all the help of classical regularity is needed. Once the appropriate, 
purely deterministic sensitivity estimates are established, it is the principle of concentration 
of measure that provides the stochastic estimates on the random solution itself. In a work in 
preparation with Hendrik Weber, the last author is applying this philosophy to the fully non¬ 
linear parabolic equation £> + dt4> — 5Jtt( 0) = £ forced by space-time white noise £ to establish 
Holder-^— bounds with exponential moments in probability. 

2 Precise setting and motivation for this work 

While the contribution of this note is purely deterministic, and the main result will be stated 
without reference to probabilities in the next section, the motivation is probabilistic and will 
be given now. In elliptic homogenization, one is interested in uniformly elliptic coefficient fields 
a in d-dimensional space W l , where uniform ellipticity means that there exists a (once for all 
fixed) constant A > 0 such that 

\/x £ M d ,£ £ £ • a(x)£ > A|£| 2 , |a(x)£| < |£|. (1) 

To such a coefficient field a we associate an elliptic operator in divergence form —V • aV. For 
simplicity (in order to avoid dealing with the correctors of the dual equation), we shall assume 
that a(x), and thus the corresponding operator, is symmetric. We note that statements and 
proofs remain valid in the case of systems with the above strong ellipticity property, but for 
simplicity, we shall stick to scalar notation like in (jT]), where we wrote £ £ W l instead of a 
tensor-valued object. 

In stochastic homogenization, one is interested in ensembles of uniformly elliptic coefficient 
fields, that is, probability measures on this space. We denote by (•) the corresponding expecta¬ 
tion and use the same symbol to refer to the ensemble. Minimal requirements for homogenization 
are stationarity and ergodicity, where both notions refer to the action of the translation group 
W l on the space of uniformly elliptic coefficient fields. Stationarity means that the distribution 
of the random field a is invariant under shifts a(z + •) for any shift vector z £ R d . Ergod¬ 
icity means that shift-invariant random variables, that is, functionals a H y £(a) that satisfy 
£ (a(z + •)) = £ (a) for all shifts z, must be (almost surely) constant. Under these assumptions, 
the classical theory of (qualitative) stochastic homogenization introduced by Kozlov [17J and by 
Papanicolaou & Varadhan m establishes the (almost sure) existence of sublinear correctors. 
More precisely, for any coordinate direction i = 1 ,... ,d and a given realization a of the coef¬ 
ficient fields, the corrector (f>i = £>j(a) modifies the affine coordinate x H > Xi to an a-harmonic 
coordinate x >Xi + 4>i{x), that is, 


- V • o(ej + V(j>i) = 0. 


(2) 


In order to be rightfully named a corrector, (f>i should be dominated by Xi, that is, have a 
sublinear growth, at least in the L 2 -averaged sense of 


lim — 

Rfoo R 


/ (*-/ - 

Jb r Jb r 


4>i) 2 dx 


= 0 . 


(3) 
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Under the assumptions of stationarity and ergodicity for (•) the classical theory constructs d 
functions fa(a,x),... ,fa{a,x) such that Q and (|3|) are satisfied for (-)-a. e. a. 

The actual construction of fa for a fixed coordinate direction i = 1 ,,d proceeds as follows: In 
a first step, one constructs a random vector field fa{a,x ) € that is stationary in the sense of 
fa(a{z+-)),x) = fa{a, x+z), of expected value {fa) = e* and of finite second moments (|/j| 2 ) < oo 
and that is curl-free djfak — dkfij = 0 and satisfies the divergence condition V • (a fa) = 0. Using 
the stationarity to replace spatial derivatives by “horizontal derivatives”, the existence and 
uniqueness of these harmonic 1-forms follows from Lax-Milgram. The second step is to consider 
the random field fa{a, x ) which satisfies 'S/fa = fa — {fa) = fa — c% — and thus satisfies (J2]) — and 
is (somewhat arbitrarily) made unique by the anchoring fa {x = 0) =0 and thus generically non¬ 
stationary. One then appeals to ergodicity (von Neumann’s version combined with a maximal 
function estimate) to show with help of Poincare’s inequality that (' Vfa ) = 0 translates into Q. 
Incidentally, the homogenized coefficients are then given by = (a(ej + S/fa)) = {afa) and 
are easily shown to satisfy (in the symmetric case) the same ellipticity bounds as the original 
ones, c. f. (HI), 

£ • cihi > A|£| 2 , k£|<|£|. (4) 

We now make the point that next to the scalar potential fa, it is also natural to consider a 
vector potential, in the parlance of three-dimensional vector calculus. Indeed, the above (purely 
functional analytic) first step constructs the harmonic stationary 1-form _/). This means that 
next to the stationary closed 1-form fa , there is the stationary closed (d — l)-form q j := afa. In 
the language of a conduction model, where the tensor a is the conductivity, fa corresponds to 
the electric field whereas afa corresponds to the current density. Hence next to considering a 
(non-stationary) 0-form 4>i with dfa = fa — {fa), it is natural to also consider a (non-stationary) 
(d — 2)-form Oi with dcq = qi — {qi), where d denotes the exterior derivative. In Euclidean 
coordinates, a (d — 2)-form a is represented by a skew-symmetric tensor {(Tjk}j,k=i,...,d> that is, 
c Tjk + crkj = 0, then dcq = qi — {qi) translates into SZ-Oi = qi~{qi), where {S/ofaj := d^o^k, using 
Einstein’s convention of summation over repeated indices. In view of the definitions qi = afa, 
fi = e % + S/fa, and ahCi = (a(ej + S/fa)), this implies 

S7 ■ (j t = a{ei + Vi fa) — ahCi and cr, is skew. (5) 


The merit of this vector potential cq has been recognized in the case of periodic homogenization 
and lies in a good representation of the homogenization error: When comparing a solution of 
—V • aS/u = f to a solution of the corresponding homogenized problem —V • ahS/v = /, more 
precisely to the “corrected” solution v + fadiV (with summation convention), one obtains for 
the homogenization error w := u — (v + fadiv) the following simple equation: —V • aS/w = 
V • {{faa — Oi)S7div). It is a that allows to bring the r. h. s. into divergence form, which makes 
simple energy estimates possible. 

In [9], we follow the classical arguments to show existence of a sublinear a under the mere 
assumptions of stationarity and ergodicity. More precisely, given a coordinate direction i = 
1 ,... ,d, we show by a suitable choice of gauge that there exists a skew symmetric tensor field 
Oi{a,x) such that its gradient Vcrj is stationary, of mean zero, of finite second moments and 
such that S7 ■ Oi = qi — {q^. We then appeal to the same arguments as for fa to conclude 


lim — 

-Rfoo I! 



( 6 ) 


for almost every realization a of the ensemble (•). 
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The contribution of this note is purely deterministic in the following sense: We will consider a 
fixed coefficient field a which is uniformly elliptic in the sense of (HD and a constant coefficient ah, 
also elliptic in the sense of ©. We then assume that for every coordinate direction i = 1,... ,d, 
there exists a scalar field <f>i with © and a skew-symmetric tensor field cr* with ([5]). Our main 
assumption will be a quantification of the sublinear growth © and ([6]), roughly in form of 

ix Y =o (7) 

for some exponent a > 0 and for (-)-almost every a. Obviously, this quantification of the 
sublinear growth of fields that have a stationary gradient of vanishing expectation can only 
be true for a quantified ergodicity of (•). The discussion and proof of criteria under which 
assumption on (•) the strengthened sublinearity (0 holds is not part of this note, but of ongoing 
work. We just mention here that for d > 2 and under the assumption that (•) satisfies a spectral 
gap estimate with respect to Glauber dynamics, it is shown in [13] , extending the arguments of 
im from the case of a discrete medium to a continuum medium, that a stationary corrector (pi 
exists, which in particular implies © for the ^-part for every a < 1. We also mention that © 
is expected to hold for some a > 0 under much weaker assumption on (•), as is already suggested 
by PI, and more explicitly in [TJ, and will be formulated as needed here in an updated version 
of [9j, where the rate of decay of correlations of a will be related to a. 

3 Main results 

In this section, we state and comment our main results, namely Theorem [2] and Corollary [3l As 
mentioned in the introduction, our standing assumption is the uniform ellipticity of the fixed 
variable coefficient a on W l and the constant coefficient ah, c. f. © and ©, respectively, and 
the existence of the scalar and vector potentials ( cp,a ) = {((pi,&i)}i=i } d } with © & ©. 

Our result, Theorem [2] and Corollary [3l on the homogenization error relies and expands on a 
large-scale regularity result established in [9] and stated in Theorem [TJ It falls into the realm of 
Campanato iteration for the C fl,Q! -Schauder theory for elliptic equations in divergence form. In 
the framework of this theory, the C' 1, “-H61der-semi-norm is expressed in terms of a Campanato 
norm that monitors the decay of the (spatially averaged) energy distance to affine functions 
as the radius becomes small. In Theorem |T] below, affine functions in the Euclidean sense 
are replaced by a-linear functions, the space of which is spanned by the constants and the d 
correctors (p±,... ,(pd- This intrinsic “excess decay” (in the parlance of De Giorgi’s regularity 
theory for minimal surfaces) kicks in on scales where the corrector pair (cp, a) is sufficiently 
sublinear, c. f. ©. This is not surprising since in view of © & ©, which we rewrite as 
(a — a/j)ej = V • cr* — aVcpi, (<p, a) is an averaged measure of distance of the variable coefficient 
a to the constant coefficient ah, whose harmonic functions of course feature excess decay. In 
this sense, Theorem [T| is a perturbation result around constant coefficients. In the context of 
stochastic homogenization, Theorem |T| amounts to a large-scale regularity result for a-harmonic 
functions, since in view of © and ©, the smallness © is expected to kick in above a certain 
length scale r* < 00 . 

The idea of perturbing around the homogenized coefficient ah in a Campanato-type iteration to 
obtain a large-scale regularity theory (Schauder and eventually Calderon-Zygmund) is due to 
Avellaneda & Lin, who carried this program out in the periodic case [3]. Recently, Armstrong & 
Smart [2] showed that this philosophy extends to the random case, a major insight in stochastic 
homogenization. This in turn inspired [9] to make an even closer connection by introducing the 
intrinsic excess and by using the vector potential a to establish its decay. 
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Theorem 1 (Gloria, Neukamm, 0.). Let a € (0,1). Then there exists a constant C = C(d, A, a) 
with the following property: Suppose that the corrector has only a mild linear growth in the point 
x = 0 in the sense that there exists a radius r* such that 

(■f I for r > r*. (8) 

\JB r JB r J L/ 

Then for any radii r* < r < R and every a-harmonic function u in {\x\ < R} we have 

i i 

inf (1 |V«-C i ( ei + V0 i )| 2 V <c(^) a inf (-f |Vn-^( ei + V^)| 2 V- ( 9 ) 

£eR d \JB r J ?£R d \Jb r J 


Here conies the first main result on the homogenization error. It relates the homogenization 
error to the amount of sublinear growth of the corrector couple (</>, a). 

Theorem 2. Suppose that the corrector grows sublinearly in two points x G {0, xo} in the sense 
that there exists an exponent a € (0,1) and there exists a radius r* < oo such that 




for r > r*. 


( 10 ) 


For a square integrable vector field g we compare Vrt defined through 


— V • aS7u = V • g (11) 

with div(ei + X7(j)i) (using Einstein’s summation convention of summation over repeated indices), 
where Vu is defined through 


- V • a h Vv = V • ( gi(ei + d ^)). (12) 

Provided suppg C B rt ( 0) and |xo| > 4r* we have 

|v "“ i>i '’ (ei+v * )|2 ) 1 ^ o ( d ’ x ’ a ) - j 0 r °{ I ls|2 ) 5 ' (13) 

Theorem [2] compares the Helmholtz projection T based on a with the Helmholtz projection T^ 
based on ah and the multiplication operator M := id + TLf. Loosely speaking, it states that 
T w M T h M*. 

We post-process Theorem [2] to get a measure of the homogenization error on the level of the 
Green’s functions in Corollary [31 our second main result. More precisely, we compare the 
“quenched” Green’s function G(x, y) and the Green’s function (or rather fundamental solution) 
Gh(x) belonging to the homogenized coefficient ah characterized by 

— V x • a(x)V x G(x,y) = 6(x — y) and -V-a h \7G h = 5. (14) 

Let us make two remarks on the existence of the quenched Green’s function: 1) In case of 
d = 2, the definition of G is at best ambiguous. However Corollary [3] only involves gradients 
of the Green’s function which are unambiguously defined and for instance constructed via ap¬ 
proximation through a massive term or through Dirichlet boundary conditions. 2) DeGiorgi’s 
counterexample [7j implies that in the system’s case, there are uniformly elliptic coefficient fields 
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that do not admit a Green’s function. However, in [5] we show that under the mere assumption 
of stationarity of an ensemble (•) of coefficients, the Green’s function G(a, x, y ) exists for (-)-a. 
e. a. We thus will not worry about existence in this note. 

The corollary compares G(x,y ) to Gh(x — y) on the level of the mixed second derivatives 
\7 x \7 y G(x,y) (interpreted as a 1-1 tensor) and —V 2 G/ l (x — y), where the mixed derivative 
of the homogenized Green’s function is corrected in both variables, leading to the expression 
—didjGh(x — y)(ei + Vfifix)) (g> (ej + V(fj(y)). The corollary monitors the rate of decay of this 
difference in an almost pointwise way, just locally averaged over x ~ xq and y ~ 0 and shows 
that, up to a logarithm, the rate of the decay is |a?o|~ d_Q , which is by |xo| _Q stronger than the 
rate of decay of |xo|~ rf of the —\7 2 Gh(x o). The main insight is thus that this relative error of 
|x 0 | _Q is dominated by the sublinear growth rate of the corrector couple (0, <r), where it is only 
necessary to control that growth at the two points of interest, that is, 0 and xq. In this sense, 
Corollary [3] expresses a local one-to-one correspondence between the sublinear growth of the 
corrector and the homogenization error. 

Corollary 3. Suppose for some exponent a £ (0,1) and some radius r* the corrector couple 
(4>,cr) satisfies 0771) in two points x = 0,rro of distance |xo| > 4r*. Then we have (using 
summation convention) 


Gn) ^B 


Go) (0) 


I V x V y G(x,y) + didjG h (x - y)(e, + Vfifix)) <g> (ej + Vfij(y ))| 2 dydx 

In^l 

< C(d, A, a) 


^ Go I 'jdA-a 


In order to pass from Theorem [2] to Corollary [3l we need the following statement on families 
(rather ensembles) of a-harmonic functions, which is of independent interest and motivated by 

Ill- 


Lemma 4. For some radius R, we consider an ensemble (■) (unrelated to the one coming from 
homogenization) of a-harmonic functions u in {|x| < 2R}. Then we have 


'M<# 


|Vu| 2 dx 


< sup(|Fu| 2 ), 
F 


where the supremum runs over all linear functionals F bounded in the sense of 


(15) 


|T\t| 2 < [ |Vu| 2 dx, (16) 

J\x\<R 

and where here and in the proof < means < C with a generic C = C(d, A). 

We note that the statement of Lemma [4] is trivial for an ensemble (•) supported on a single 
function. Hence Lemma |7] expresses some compactness of ensembles of a-harmonic functions. 


4 Proofs 

Proof of Theorem [2J Throughout the proof < denotes < C , where C is a generic constant 
that only depends on the dimension d, the ellipticity ratio A > 0, and the exponent a £ (0,1). 
By a rescaling of space we may assume w. 1. o. g. that r* = 1; by homogeneity, we may w. 1. o. 
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g. assume (f \g\ 2 ) 2 = 1. We set for abbreviation R := \\xq\ > 1. We first list and motivate the 
main steps in the proof. 

Step 1. In the first step, we upgrade the excess decay Q in the sense that we replace the 
optimal slope £ on the 1. h. s. by a fixed slope that does not depend on r, but is the optimal 
slope for some radius ro of order one. More precisely, suppose that (fTUl) holds for a point x , say 
x = 0. Then there exists a radius 1 < ro < 1 such that for all radii R > r o and a-harmonic 
functions u in {|x| < R}, the optimal slope on scale ro 


£ := argrnin 



I + V<(>i )| 2 


is such that for r G [ro, R]: 

i i 

{L |Vn ■ Uei+ v ^ )|2 ) 2 ~ { i )i+a ( L |vu|2 ) 2 ■ (17) 

In addition, we have 

l£l < (£ |Vn| 2 ) 2 < (J b |Vu| 2 ) 2 . (18) 

Since rg < 1, the second inequality holds, possibly with a worse constant, for all 1 < r < R. 


Step 2. We have the following estimates on Vu and Vu 

(/l v “| 2 ) ;$1, sup| x |> 2 (|x| d |Vu(x)| + |x| d+1 |V 2 u(x)|) < 1. (19) 

Moreover, u and v have vanishing “constant invariant” 


V? 7 r • aS7u = / V?? r • ah'Vv = 0 


( 20 ) 


and, thanks to the special form of the r. h. s. of (1121) . identical “linear invariants” (for k = 

j Vr/ r • {{x k + (t> k )dVu - ua(e k + V(j> k )) = j Vr/ r • (x k a h Vv - va h e k ), (21) 

for all r > 1, and where rj r (x) = rj(^) and r/ is a cut-off function for {|x| < 1} in {Jx' < 2}. 
We speak of invariants, since for two a-harmonic functions u and u (in our case u = 1 for the 
constant invariant and u = x k + (j) k for the linear invariant) defined in {|x| > 1}, the value of 
the boundary integral f. )(} v ■ ( uaVu — uaVu ) does not depend on the open set Q provided the 
latter contains B\. 


Step 3. We consider the homogenization error 

w := u - (v + cfridiv), 

where cj) denotes the following blended version of the corrector: 


< i-= (l-vM- 


B i(0) 


4) + v(<P- T 

Je 


B i(a-o) 


0), 


( 22 ) 


(23) 


where is a cut-off function for {|x — xq\ < i?} in {\x — xq| < 2 R}. Then we have 


— V • a\7w = V • h for \x\ > 2, 


(24) 
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where Vtc and h satisfy the following estimates 


'M>2 


|Vw| 2 <1 


(25) 


and 


' |a;|>r 


\h \ 2 <~ 


T 2 


+a 


for r > 2, 


\h\ 


2 — 

2 1 < r2 
~ J^d+c 


for 1 < r < R. 


V J \x—xq\<t 

Moreover, w has an asymptotically vanishing constant invariant 

lim / Vr] r ■ a\7w = 0 

rtoo J 

and asymptotically vanishing linear invariants 

lim / \7r) r ■ ((xk + ^)aVio — wa{ek + Vr/)^)) = 0 for k = 1,..., d. 
rtoo J 


(26) 

(27) 

(28) 

(29) 


Step 4. Extension into the origin. There exists w (with square integrable gradient), h, and / 
defined on all of M. d such that 

— V • aS7w = V • h + /, (30) 

while 

Vw = Vw in {|x| > 4}, supp h C {|x| > 2}, supp / C {|x| < 4}. (31) 

The r. h. s. / and h satisfy the estimates 


P) <1, 


(32) 


and 




\hf 


< _ 

~ ± +a 
j, 2 ' ^ 


I JP 


d 

r 2 

< — _ 

~ Rd+c 


for r > 2, 


for 1 < r < R. 


\J \x—xo\<r 

Moreover, we have vanishing constant and linear invariants: 

Jf = 0, J ((e fc + V0 fc ) • h - (x k + 4>k)f) = 0 for k = 1,... ,d. 

The estimates p2]) - (lMll ensure that the integrals in (1351) converge absolutely. 

Step 5. From the equation (1301) . the vanishing invariants (1351) . and the estimates 
it follows 

i 

Ini? 


(33) 

(34) 

(35) 


& 


-12 


\x\>R 


|V«)| 


< 


i?2 


+Q 


(36) 
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Step 6. From the same ingredients as in Step 5 and in addition (1341) . one obtains the following 
localized version of (l36l) : 



< 


In R 
R d+a ' 


Step 7. Conclusion. 

Argument for Step 1 

We first argue that there exists a radius 1 < ?’o ~ 1 such that 


Ifi ' (e* + V ^)| 2 ~ |f | 2 for r > r 0 . 


B r 


(37) 


Indeed, (1371) easily follows from the sublinear growth (fTOl) of p in the point 0. The upper bound 
is a consequence of Caccioppoli’s estimate for the a-harmonic function Xi + (pi — f Bo pi), cf. 
(P|). leading to 



|& ' (e* + V0j)| 


1 

2 



I £i( x i + 




The lower bound is a consequence of Jensen’s inequality once we introduce a cut-off function r] 

of B i in B\ and set r] r (x) = 77 (f): 

2 r 



If* • ( e i + 



> / rj r \ii ■ (ei + Vpi 


> 


1 


(fVr)? 

1 

(fnr)i 
C 


J Vr£i ■ (ei + VPi) 

f Vr - fi (Pi ~ T Pi)Vri r 

J J J Br 


> If I ( 7^2 -CA 1 ( ( \Pi- i 
\Jb t Jb 


B r 


For an a-harmonic function u in Br and a radius tq < r < R we consider 



fr := argmin -r Vu - fj(e; + Vpi) 2 

J B r 

(38) 

and claim that 

1 



If, - 6 ?| < (inf £ |Vu - Uei + Vpi) | 2 ) ' • 

(39) 


Indeed, by the triangle inequality in M. d and since a > 0, it is enough to show for r$ < r < r' < R 
with r' < 2 r that 


Ifr - fr'l < (^)“ (inf £ | Vu - ti(ei + V&) | 2 ) * , 
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which thanks to Theorem [T) follows from 

l^-Cr'l 2 <mf/ \Vu-^i(ei + V(f>i)\ 2 + m£-f |V« - &(«* + V0;)| 2 , 

? £ J B r i 

which by definition (1381) in turn follows from 

|£r-£r'| 2 < / |Vri-?r-,j(ei + V(/>i )| 2 + / |V« - £r',i( e i + V<^i)| 2 . 

The latter finally follows from the triangle inequality in L 2 (B r ) and from (1371) applied to £ = 

Now (fT71) is an easy consequence of (f39l) (with (r,R) replaced by (ro,r)) and Q: 

/ r \5 1391 ( r _ \ 2 19} r , / r \ \ 

(J b \Vu-ae l + V^) < (inf |Vu-&(e* + V0*)| 2 j < (-)^+ Q ^ |Vu| 2 J . 


We finally turn to the argument for (1181) . For this purpose we first note that for all radii r > ro 


16-1+( mf £ |Vu-^( ei + V^)| 2 


' B r 


< 

n-j 


/ |Vu|‘ 

J jBf 


(40) 


Indeed, that the second 1. h. s. term is dominated by the r. h. s. is obvious. For the first 1. h. s. 
term we note by (1371) . the triangle inequality in L 2 (B r ), and the definition (1381) of that 


l&l <(/ l^(e* + V^)| 2 V < (inf/ |Vu-^( ei + V^)| 2 V + f/ |Vu| 
\4B r / V ? JB r J \JB r 


Equipped with (ITUl) . and more importantly (1391) and ©, we may now tackle (fl 8 l) . For the hrst 
estimate in (fT 8 l) . we appeal to the triangle inequality, (1391) with ( r,R ) replaced by (ro,?’), and 

(®D: 


l£rol — |£r| T |Cro £r \ 


ICrl+finf / | Vu-^ + V^pY ^ (/ |V«| 
V ? J B r J \JBr 


For the second estimate in (1181) . we use the triangle inequality in L 2 (B r ) and in K d , 
and (HOI) : 


(f Br l V-u l 2 ) 2 - l&*l + \€r - Crl + (inff f Br |Vu - ii{ei + V^)| 


< \£r\ + max p=r,R [ inf € f Bp |Vu - &(ei + Vfa )|‘ 

® / P 

< | 6 ?| + (inf ? f B \Vu - &(«* + V0*)l 

SI, 

S (h„ iv»l 

Argument for Step 2. The first estimate in (j!9l) is an immediate consequence of the definition 
dill) and the energy estimate. We now turn to the second estimate in (1191) . We first note that 
the r. h. s. g := f/,;(e* + di4 >) of (fT2l) can be bounded with help of Caccioppoli’s estimate 


\9\ < / 9: 




I ei + V0* 


< 

nsj 


\x + (6 - 


'b 2 


( 0 - / 0)1 

Jb 2 


< 1 . 
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Since moreover, it is supported in B\, we obtain the following representation for \x\ > 2 

Vv(x) = J V 2 G h (x - y)g(y) d y, V 2 v(x) = j V 3 G h (x - y)g{y) d y 

in terms of the constant-coefficient Green’s function Gh■ The second estimate in m now 
follows using the homogeneity of Gh- 

We now turn to the invariants. The first identity (1201) follows immediately from integration 
by parts of the equations (fill) k (fl2l) using the fact that the respective right hand sides are 
supported in B\. We now turn to the second identity m, which also follows from integration 
by parts, but this time relying on the Green formulas 

V • [{x k + fa){aVu + g) - ua(e k + Vfa)} = (e k + Vfa) • g, 

V • [x k {a h Vv + gi(ei + di<j>)) - va h e k ] = e k • [^(e* + dfa)], 

which follow from ED k (fT2l) in conjunction with and the pointwise identity 

(e k + Vfa) ■ g = e k • fafa + 8$)]. 

Argument for Step 3. We start by establishing the formula (fMl) with 

h := (4>ia - a^VdiV + div((f fa ~ T fa) a - (f f cr») > ) V 77 , (41) 

V 4Bi(0) JB 1 (x 0 ) JBi_( 0) JB^xq) J 

where, in line with (1231) . we have set 

a := (1 — rj)(a — -r a) + g(a - -r a). (42) 

J Bi{0) J Bi(xo) 

Indeed, from definition (1221) k (1231) we obtain 

Vw = Vu - [divfai + Vfa + (f fa-f fa)V g) + faV dfa (43) 

J Bi(0) J Bi(xo) 


and thus, using V • afa + Vfa) = 0, cf. (j2j). 


—V • aVw = —V • aVu + V<9 iv ■ afa + Vfa) + V • [div(-r fa — -jr fa)aVg + faaVdj,v\. 

J Bi(0) J B\(xq) 

Using the identity afa + Vfa) = 4- V • cq, cf. ©, we have Vdtv ■ afa + Vfa) = V ■ a,hVv + 

VdiV ■ (V • ai). Using that the r. h. s. of (TTTT) k (fT21) are supported in B 1 , we thus obtain in 

{\x\ > 1 } 

— V • aVw = VdiV ■ (V • <7j) + V ■ [e?iu(-r fa — -r fa)aVi 7 + faaVdiv]. 

J Bi(0) J B\(xo) 

It remains to substitute a by a and to bring the related part of the above r. h. s. into divergence 
form. Indeed, by definition (|42l) 

V • a* = V • di - V ■ (g(f ai-f <Ji)) = V ■ fa - (f °i-T <?i)Vri, 

J B\(0) J B\(xq) J Bi(0) J Bi(xo) 

and by the identities 


VC ■ (V • a) = V ■ (CV • a) = -V • (trVC) for skew a, 


(44) 
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we obtain as desired 


Vdiv ■ (V • o-j) = -V • [diVdiv] - V ■ [div(f cq - / 

J .Bi(O) JB-\ fxn) 


' Bi(xq) 


In order to estimate the contribution to Vw and h that comes from the difference f Bl (o) — f_B 1 (x 0 ) 
of the average values we now argue that 


/ (<^0-) -/ (0,cr) 

JB i(0) J Bi(xn) 


B 1(10) 


<R 


1 —a 


(45) 


To keep notation light, we write </> instead of (</>, a). Let us first argue how to reduce (1451) to 


/ *"/ 

J B r (x) JB 


Bx(x) 


< r 1 - 01 


for r > 1, x £ {0, xq}- 


(46) 


Indeed, (|45l) follows from (PT6l) and our assumption (fTOl) via the string of inequalities 


/ 0 “/ 

J -Bi(O) JBi Oo) 

< (/ x 10-/ 0 i 2 V +(7 , 10-/ 01 s 

\Jb 2R (±xo) Jb^O) J \Jb 2R (\x 0 ) iBi ( xo ) 

< (/ 10 -/ 0l 2 V + f/ 10-/ 0l 2 V 

\JBar( 0) JSlfO) / \JBar(x n ) JBi ( xn ) 


< 


/ 10 -/ 0i 2 V+(7 10 -/ 

J Bar(0) J Bar(O) / \JBar(x n ) JE 


' Ba R (Q) 


Bar(xq) 


+ 


/ 0 -/ 0 + / 0 -/ 

J B4,r(0) JB 1 CO ) J Bar(x 0 ) JE 


'Bi{ 0) 


0| 2 

Bar{xq) , 

(flOl). (I46i) 

< 


' B iR (xo) 


' B 1(10) 


1—a 


For AMD , we focus on x = 0 and note that by a decomposition into dyadic radii, it is enough to 
show 


/ 0-/ 

J B r J -B_/ 


< r 1 a for 2r > / > r > 1. 

This estimate follows from a similar string of inequalities as the one before: 


/ 0-/ 

J B r J £L/ 


< 


< 

r^j 


/ 10 - / 01 " 

J JB'p J ^B 'p 

/ 10-/ </| 

J jBr J jBr 


+ 


+ 


/ 10 -/ 0| 2 

^ J B r J B r f 

/ 10 -/ 0| 5 

, JB, JB. 


and an application of (fTOl) . 

We now turn to (|25l) . We start from the formula (1431) . which we rewrite as 


\7w = Vu - diviei + Vfa + (t 0; - 7 0j)Vi?) 

V J Si( 0 ) JBi{xo) 7 

-/ 0O + 7 ?(/ 0* — / (j>i))vdiV. 

JBi(O) JBU 0) iBihn) 7 
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( 47 ) 


From the estimate (I45|) on averages and the estimate (1191) on v we thus obtain for |x| > 2 
\Vw\ < |Vn| + |x| _d (|id + V<f>\ + i? 1-Q |Vr/|) + |d _(d+1) (\<f> - / <f>\ +R 1 ~ a r]). 

JBRO) 

Taking the L 2 ({|x| > 2})-norm, we see that the V'«-term is bounded according to (1191) . By the 
choice of the cut-off function 77 , the function |a:| _rf |V 7 /|-|-|x |^ ( '' i+1 ^77 is supported in {|x—xo| < 27?} 

and bounded by so that the contribution of this term to 


|x | >2 l Vw l J is estimated 

by R~ a ~ 2 . We turn to the term involving 0 and will for later purposes show the slightly more 
general statement for r > 1 


' |:r|>r 


(1*1 (d+1) l(0w) - / (<t>,°)\) 

J Bi(0) 


1 


~ A+ a ' 


(48) 


which, restricting to <f> in our notation, follows by dyadic summation from 


> r<\x\<2r 


10 -/ 01 

Jbro) 


2 < _§+l-a 


which trivially follows from 


1 

/ 10-/ 0| 2 )<r 1 - a 

JBr( 0) JBR 0) / 


The last estimate is a combination of our assumption (1101) with (I46|) . We now turn to the 
estimate of the V0-term in (1471) . By dyadic summation, it is enough to show 


'L 


I id + V0 | 2 <r 2 , 


I r<\x\<2r 

which follows from Caccioppoli’s estimate and assumption (HOD: 

(/ lid + V^I 2 ') <-(/ | x + (0-/ 0)| 2 ^) <! + -(/ 10-/ 0| 2 ) 

\J i?2r / i?4r 7 £?4r / ^ \7 £?4 r J B^ r J 

We now turn to estimate (1261) . For this purpose, we rewrite the definition (1411) of h as 


h = 


l 


a - <7,; - 


B i(0) 

4*i 

B\(0) J Bi(xq) 


£ 


B i(0) 

/ 0*) a — (/ cr* - T CTi)) 

JSifxo) 7 Bi(0) IBi(id) 7 


a,-) Va-v 


+ ^( 

+ d i v(('f 0 ; - / 4 >i)a — (/ 

v 7 Bi( 0 ) JBi(xn) 


o’*) IV dj,v 


B i(0) J B\(xo) 


&i - 
Bi(0) J Bi(x o) 


/ ch))v??. 

7Si(xo) 7 


Inserting, as above, the estimate m on v and the estimate (|45D on the averages we obtain for 

Id > 2 


\h\ < |x|- (d+1) ( 0 ,( 7 ) - / ( 0 ,< 7 ) +R 1 - a (\x\-^7 1 + \x\- d \\/r ] \). 


B i(0) 


(49) 
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Since by definition of rj, the function |x|~( rf ++ + |u;| —rf |V 77 I is supported in {|x| > 2R} and 

bounded by R~( d+1 \ its contribution to |/i| 2 ^ vanishes for r < R and is bounded by 

R 2 R-( d + a ) < r ~(|+ a ) for r > R. The first r. h. s. term in (|49l) was treated in (l48l) . 

We finally turn to the last estimate m- For this purpose, we note that because of the choice 
of r/, on Br(x 0 ), the definition (TTTI) of h turns into 


h = 


BRx 0 ) 


4>i)a - (ai - 4- (Ti))VdiV, 
Jb^xq) ' 


so that by the estimate m on v we have for 1 < r < R 


/ \h\ 

B r ( xq ) 


2 j < R~( d + 1 ) 


/ I ( 0+1 

J B r (x 0 ) J Bi(x 0 ) 


which by estimate (14611 on averages turns into 

(i I (0+-/ (0+l 2 l + + 

\ J Br(xn) / \ \JB r (x 0 ) J B r (x 0 ) ) 


so that the desired estimate in the strengthened form of 

1 

( r \ 2 

\2 \ 


+ 1 — Q! 


I B r ( xo) 


I hY 


r 2 

< - _ 

~ R d+1 


now follows from assumption ( 1101 ) . 

We finally turn to the asymptotic invariants (1281) &: m- Since we may assume r > 6 R, we may 
ignore the presence of r] in the definition (1221) of w (and assume w. 1 . o. g. that f Bi ^((j),a) = 0 
for notational simplicity). Hence we have 


aS7w = aVu — ( dj,va(ei + V^i) + 4>id\7diV ) ^ aVu — (o^Vu + c+V • < 7 * + fadS/div). (50) 


Using the formula 


V?y r • (CV • a + erVC) = 0 for skew a, 


(51) 


which follows from the last identity in (14411 (and even holds if £ a priori is not defined in {|x| < r} 
since it can be arbitrarily extended), we derive 


J Vr) r ■ aVw = J \7rj r ■ (aS7u — a^Vr — + a — <7i)Vc+), 

which by the identity ( 1201 ) of the variants for u and v collapses into 

J \7r] r ■ aVw = — J V?? r • (+0 — ai)X7di.v ). 

Using the estimates (fTUI) on v, this yields 


V?y r • a\7w 


< r ~ 0+2) 


/ 

J B2r 


+)l ^ r 


-2 


/ l(0+f 

J B2r 


(52) 
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Together with our assumption (fTU|) . this implies (I28f) . 

We conclude this step with the argument for (1291) . From the identity (1501) we deduce 


(x k +cj) k )aVw - wa(e k + V^fc) 

= ((x k + (/>fc)aVn - ua(e k + Vfa)) - ( x k a h v - va h e k ) - (j) k a h \7v 
- (x k + 4> k )(divV ■ Oi + faaVdiv) + <j>idiva(e k + V0 fc ) + nV • a k . 

Using that the linear invariants of u and v coincide, cf. m, and the formula m, only the 
following terms survive after application of fVr/ r -: 

-(frkCLhVv + o-,;V((.Tfc + (t>k)div) - (, x k + dw + (j)idiva(e k + V^) - cr k Vv 

= div(4>ia + ai)(e k + V0 fc ) - (4> k ah + er fc )Vu - (x k + <A0(<M - o‘i)V5 i u. 

Together with the estimates (fT9l) on v, this implies by Caccioppoli’s estimate 
/ Vr] r ■ ((xk + 4>k)a\7w - wa(e k + V0 fc )) 


<r (d+1) f | 0 ,cr)|(l + |id + W<f>\) + r (d+2) f a)\\x + <j>\ 
J B“2r J B“2r 



By assumption (fTOl) . this yields (l29h . 

Argument for Step 4. Select a cut-off function 77 for {|x| > 4} in {|x| > 2} and set 
w := rj(w — wq), h := r/h — (w — Wo)ciVr), f := — V 77 • (a Via + h), 


where wq is the average of w on the annulus {2 < \x\ < 4}. By the choice of 77 , (13111 is clearly 
satisfied. Since 


Vw = r/V'w + (w — wo)Vr) and thus 


— aVw — h = —77 (aVzu + h), 


we learn from (1241) that also (1301) holds. The estimate (1321) on / follows from the estimates (1251) 
and (1261) on w and h. As for estimate (1331) on h, we note that for the contribution 77/1 to h, the 
estimate immediately translates from (|26l) : for the contribution — {w — t«o)aV) 7 we note that it 
is supported in {2 < \x\ < 4} and estimated by \w — wo\, and thus the desired estimate follows 
from Poincare’s inequality with mean value zero on the annulus {2 < \x\ < 4}. For estimate 
(1341) on h, we note that on Br(x 0 ), h coincides with h, so that it follows immediately from (1271) . 

We now turn to the invariants (1351) . For the constant invariant we note that we obtain from the 
equation (1301) for r > 2 




, _ j. 1311 

Vr/ r ■ (.aVw + h) = 


/ 


\/7] r ■ (aVw + h). 


Hence the first identity in (1351) follows from (1281) provided we have lim^oo J \7r] r ■ h = 0. The 
latter is a consequence of ( 1261 ) : 
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We now turn to the linear invariants in (1351) . We start with the identity which follows from the 
equations for w and 4>k and the support properties 


/ 


Vr( ~ (efe + V^fc) -h + (x k + (t>k)J) 

J30I f 

= / V(?/ r (xfc + (j>k)) ■ a\7w + (x k + ■ h 

j2j r 

= / Vry r • ((x fe + (j) k )(aVw + h) - (w - w 0 )a(e k + V^*)) 
= / V?7 r • ((xfc + cj)k)(a\7w + h) - tca(e fc + V</>*,)) • 

Hence the linear invariants follow from (1291) once we show 


- (e k + V^fc) • h + (x k + 4>k)f\ < oo 


and 


lim / (x k + cP k )Vrir • h = 0. 

rtoo ,/ 

The limit (1541) follows from the estimate 


(53) 

(54) 


(Xk + 4>k)^1]r ' h 


< \A + - ( [ | A 12 

\x\<2r 


' |x|>r 


— 11 + 

r a \ r \./ B2r 


H /, 


together with the observation that the term 


r ( $B 2r 


(55) 


is of higher order since w. 1. o. g. 


we may assume f B <p k = 0 so that in view of (1461) we may appeal to our assumption (1101) . 

We now turn to (1531) : in view of the square integrability of / and h (established above) and the 
local square integrability of Xk + 4> k and its gradient, it remains to show 


/ \(e k + V<4) ■ h\ < oo. 

J\x\> 2 

To this purpose, we divide into dyadic annuli and use Caccioppoli’s estimate 


(56) 


I \(e k + V4>k) ■ h\ 

J r<\x\<2r 


<1/ |e fc + V0 fc | 2 / \h \ 2 

B 2r J |x|>r . 


< 

r^j 


[ Hkf 

J JDAr 


’ |a;|>r 


\hf 


. 1 / d 

^ - k 2 + 

r \ \.ib 4 . 


' |a;|>r 


I hf 


We now appeal to the same argument as for (f55l) to see 


[ |(efc + V0fe) • h\ 

J r<\x\<2r 


< 

r\_/ 
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Summation over dyadic annuli yields 

Argument for Step 5. We give an argument by duality and therefore consider for arbitrary 
square-integrable vector field h supported in {|x| > R} the finite energy solution of 

— V ■ aVw = V • h. 


Since w is a finite energy solution of (13011 we have 


/ xw ' / 


— J Vw ■ h = — I Vw ■ aVw = / (Vw ■ h — wf). 

Recall Step 1 and consider 

£ := argrnin / | Vw - £fc(e fc + V(f> k )\ 2 and c \= jr (w - £ k (x k + <t> k ))- 

J Bj'q j -B4 

By the vanishing invariants (1351) we may post-process (1581) to 


— / Vw • h = 


Vw - i k {e k + V4>k)) -h-{w- + (f> k ) - c)fj. 


By the support conditions (f5TI) on / and h, this implies 
I Vw ■ h 


< 


OO r. 

vi 

n= 1 J2 


=1 J 2 n <\x\<2 n + 1 


|Vu> - £ fc (e fc + V4> k )\\h\ + \w — + <j) k ) - c\\f\ 


< 


71= 

OO 

E 


|#|<4 


_1 B 2 n +1 


\Vw -£ fe (e fe + 


'|a;|> 2 TI 




+ ( / 1^ - + <t>k) - c[ 

J B 4, 


l/T 


(57) 

(58) 

(59) 


Inserting the estimates (1321) and (1331) on / and h and using Poincare’s inequality on £>4 (in view 
of the definition (1591) of c) 


/ 


Vu • h 


OO - 

<E<^ 

71=1 


n ' 


/ | V7Z) — £fc(cfc + V^>fc)|‘ 

«/ B2 n 


(60) 


We now distinguish the cases of 2 n < R and 2 n > R. In case of 2 n < R, since by assumption h 
is supported in {|x| > R}, w is a-harmonic in Br, cf. (1571) . We thus may appeal to (1T71) in Step 
1 and obtain in view of the definition of £ in (1591) 

(J b |Vu)-£ fc (e fc + V £> fc )| 2 ) 2 <(^)“(£ |Vu)| 2 ) 2 . (61) 

In case of 2 n > R, we use (1T81) in Step 1 (and once more Caccioppoli’s estimate in conjunction 
with assumption dlO]) ) to obtain 






Using the energy inequality for (1571) in form of 




|Vt5| 


i 

2 \ 2 ^ —- 
“ 1 <7-2 

rs _/ 


|Vu)| 


1 

2\ 2 / 
z ' N f 2 
r\j 


L 


we obtain in either case 

\\7w - £, k (e k + Vcj> k )\ 2 
Inserting this into (1601) we obtain 

Vw • h 


(£)“(£)* for 2 n < R 

(i) f 


for 2 n > R 




I 


^ Ini? 

~ R a+ i 


\h\‘ 


Since the only constraint on h was that it is supported in {|x| > R}, we obtain (1361) . 
Argument for Step 6 . In view of (1361) in Step 5, it is enough to show 


-12 


' B\{xo) 


|Vw| 


_Ini? 

<-h 

~ i?d+c 


Ml 

i?2 \Je 


-12 


Br(x q) 


|Vw)| 


(62) 


where w. 1. o. g. we may assume that i? = 2 N is dyadic. To this purpose, we decompose the r. 
h. s. h of (1301) into (h n ) n =o,...,N+ 1 such that ho is supported in B\(xo), h n is supported in the 
annulus { 2 n ~ 1 < \x — xo| < 2 n } for n = 1,..., N, and Jin+i is supported in the exterior domain 
{|x — xo\ > i?}. For n = 0,..., N let w n denote the finite energy solution of 

— V • aVw n = V • h n ] (63) 

For n = N + 1, wjy+i denotes the finite energy solution of 

- V • a\7w N+ 1 = V • h N+ 1 + /. 

By uniqueness of finite energy solutions of (1HUD . Yln=o = Vu), so that by the triangle 

inequality 


' B 1 (x 0 ) 


|Vtt>| 


2 JV+1 


n 


=0 V B i( x o) 


|Vw 


- 12 
n | 


(64) 


In the sequel, we will estimate the contributions individually. 


We start with the intermediate n = 1,... ,N: Since by construction of h n and by (1U51) . w n is 
a-harmonic in Bon-i (xo), we have by Step 1 (the second inequality in (1181) and with the origin 
replaced by xg), and the energy estimate for 


'Bi(xo) 


I Vw r 


< 


< 


' 1 \- 
—) 2 

'2 n 


1 

- ] 2 


|Vm 


- 12 


B 2 n —1 (-To) 

[ \h\ 2 

>B 2 n(x 0 ) , 
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In the case of n = 0, we obtain likewise by just the energy estimate 


|Vh) 0 | 2 < 


\J B i(xo) 

Hence in both cases we obtain thanks to (13411 that 


' B\(xo) 


\h? 


' Bi(xq) 


|VtU 


' T ’n| 2 


< - — 

~ Rd+c 


We finally turn to n = N + 1 and obtain like for n = 1,..., N by Step 1 


|Vu>AT + i|" 


^ , 1 \<L 
< ( — ) 2 


\J B\[xq) 

But now we use the triangle inequality in form of 


' B r {x o) 


| V icjv+i 1^ 


' B r {x o) 


| Vwjv+i I < 


/ |Vts| 2 V + (/ 

J B r (x o) J 


+ / |V(w - tt>jv + i)|' 


and the energy estimate for —V • aV(w — wn+i) = V ■ (I(B B (xo))h) to conclude 

1 


|V(w - twjv+i)| < 


' B r (x o) 


H 


< 


Re¬ 


collecting these estimates on the contribution of n = N + 1, we obtain 

i 

1 


' Bi (zo) 


\Vw N+l \' 


< 

~ R d+C 


+ (h~ 2 ([ l v ^l 2 


( 65 ) 


( 66 ) 


Inserting (1551) and (1551) into (1541) . we obtain 

Argument for Step 7. Note that by (f5TT) and (1451) we have in Bi(xq), where the cut-off 
satisfies r) = 1, 


X7w = Vrt - [div(ei + V<fo) + (^ - 


B 1(2:0) 


<f>i)Vdiv]. 


In 


fact, the extra term {<f) l — f Bl ( XQ \ (fr^'VdiV is of higher order: 


1 B 1 (x 0 ) 


&)V3, 

4Si(z 0 ) 


<_ 

~ Rd +1 ' 


which follows immediately from the assumption cnD and from the pointwise estimate (11911 on 


Proof of Corollary [3j Like for the proof of the Theorem [2] we may assume r* = 1. We 
consider g, u, and v like in the statement of Theorem [2] and note that we obtain from (lllH the 
Green’s function representation 


u{x) = ~ J V y G(x, y) ■ g(y) d y, 
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whereas m yields 

v(x) = j VG h (x - y) • (e* + di4>(y))gi(y) d y = j djG h (x - y)(e.j + V<^(y)) ■ g{y) d y. 
By differentiation in x this implies 

Vu(x) = — / \' x V y G(x,y)g(y)dy, 

div{x) = / didjG h (x - y)(ej + V^-(y)) • g{y) d y, 

JBi(O) 


so that (1131) takes on the form 


(Jb^xq) 



(v x V y G(x,y) + didjG h (x-y)(ei + V^(x)) 0 (ej+ V(pj(y))^g(y) 


Ay 



1 

2 


< 


ln|x 0 | 

x 0 \ d+a 


( [ \g\ 2 dy)K (67) 

JB i(0) 


We now argue that in (1H71) . we may replace didjGh(x — y) by didjGh{x). Indeed, because 
of |x 0 | >4, \x — x 0 | < 1, and \y\ < 1, we have for the constant-coefficient Green’s function 
| didjGh(x — y) — didjGh(x)\ < -p^p+r- In addition, we have by the argument from Step 1 in 
the proof of Theorem [2] and our assumption (ITOl) that f Bl ^ XQ ^ \&i + V</>j(x) \ 2 Ax < 1 as well as 
f B i(o) \ e j + | 2 Ay < 1. We therefore obtain by Cauchy-Schwarz’ inequality in y 



y) - didjG h (x )) [e.i + V^(x)) 0 (ej 


+ ^<f>j{y))g(y)dy 


2 



1 

2 



Hence (1671) upgrades to 


B 1(2:0) 



(v x V y G(x,y) + didjG h (x)(ei + V^(x)) 0 (ej + V0j 


Ay 



1 

2 


< 


ln|x 0 | 

x 0 \ d+a 


( [ \g?dy) 2 . (68) 

JB i(0) 


We now apply Lemma[3]to the family u(y) = \7 x G(x,y) — didjGh(x)(ei + Vcj>i(x))(yj + 4>j(y)) 
of Revalued maps defined for y £ -Bi(O) and parameterized by the point x £ B i(xp), i.e., we 
interchanged the role of x and y. We note that these maps are component-wise a-harmonic 
on {|y| < 1} because of (J21) and because of the y-derivative of (fill) in conjunction with the 
symmetry of the Green’s function, which follows from symmetry of A (for nonsymmetric A one 
would apply Lemma 0] for the adjoint problem). The ensemble average on this family is given 
by the spatial average (•) = / b 1 (x 0 ) "dx. The role of the linear functionals in the statement of 
Lemma [His played by Fu := f Bl / 0 ) V y u(y)g(y) Ay, where we restrict to g that are normalized 
(Jb i(o) |y| 2 dy)2 = 1. Hence we learn from Lemma [H that (1681) implies 


V 


Bl(x 0 ) 




■V y G(x,y) + didjG h (x)(ei + V(fo(x)) 0 (ej + Vcj>j(y))\ 2 dy Ax') 


< 


In |x 0 | 
\x 0 \ d+a ' 
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By the preceding argument, we may substitute didjGh(x) again by didjGh(x — y). 

Proof of Lemma [4j Using translation and scaling invariance, an elementary covering ar¬ 
gument shows that is enough to establish (fT5l) with radius replaced by -^= { . Therefore, it 

suffices to show the result with the inner ball {|x| < replaced by the cube (—^) d and 

the outer ball {|x| < R} in (fl6l) by the cube (—R,R) d . By scale invariance, we may reduce to 
(—|, j) d and (—f, §) d , respectively. We thus will show that 

([ |Vtt| 2 dx) < sup (\Fu\ 2 ), (69) 

where the supremum is taken over all linear functionals F satisfying 

\Fu\ 2 < [ |Vu| 2 dx. (70) 

k- f,|) d 


The proof essentially amounts to a generalization of Caccioppoli’s estimate for an o-harmonic 
function u. Recall that the latter states 


|Vu| 2 dx < inf 


|u — c| 2 dx, 


J(- 2 L 2 LW centj/l IW 

y 4’y 2’2> 

which we may re-express in terms of the Fourier cosine series 


Fu(k) := 



(-2L 2L\d 

v 2 ’ 2 / 


u(x) nf =1 cos (kiXi) dx for k € Z d \ {0} 


(71) 


(72) 


as f, 7T 2 L) d |Vu| 2 dx < X^fcez d \{o} \J r u(k)\ 2 . The generalization of (TTTT) we need is that for any 
even n E N we have 

|Vu| 2 dx< ]T J-\Fu(k)\ 2 , (73) 

4’V d kez d \{ o} 

which amounts to replace the L 2 -norm on the r. h. s. of (17111 by the negative H~ n - norm. For 
the remainder of the proof, < will mean up to a constant also depending on n; but this will not 
matter, since we will presently fix an n in terms of d. 

Before we give the argument for urn let us argue how to conclude. We first note that for 
k € r L d \ {0} the linear functional 

F^u := \k\Fu(k) 

has the boundedness property (I70|) . Indeed, we obtain from integration by parts in (17211 



Fu{k) 



d\u(x) sin(/ciX'i )nf =2 cos (kjXi) dx, 


so that k\\Fu{k)\ 2 < f t _n E)d (diu) 2 (x) dx. Hence, after taking the ensemble average, we may 

_ J v 2 ’ 2 > 

reformulate (1731) as 


([ |Vu| 2 dx)< 

I <_TL IW - 

y 4 * 4 J kez d \{ o} 


1 

\k\ 2 ( n+1 '> 


( \F k u \ 2 ). 


Now picking n € N with n > | — 1 so that z d \{o} | fc | 2 (i+i) < 1, we obtain (1691) . 
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We now turn to the argument for (1731) and introduce the abbreviation || • || for the T 2 ((— ^) d )- 
norm. The main ingredient is the following interpolation inequality for any function v of zero 
spatial average 


fcez d \{ 0 } 1 1 


( E i^mi 2 ) 

fcez d \{ 0} 1 1 


(74) 


where rj is a cut-off function for (—-|, \ ) d in (—f, %) d with 


IV77I < 1 . 


(75) 


Note that (174[) couples the degree of negativity of the r. h. s. norm to the degree of degeneracy 
of the cut-off rj n . If we plug the standard Caccioppoli estimate in its refined form of 

J7SJ 

\\r] n+1 Vu\\ < inf \\(u - c)Vr ? n+1 || < inf ||r f(u - c) || 

cGM cEK 

into (1741) for v = u — c and use Young’s inequality, we obtain (1731) . 

In preparation of its proof, we rewrite (1741) without Fourier transform, appealing to the repre¬ 
sentation of the Laplacian — Atv with Neumann boundary conditions through the Fourier cosine 
series by J r (—A]y)w(k) = \k\ 2 J r w(k): 

( |-7j 2^|7Y(A;)| 2 ') 2 = ||w|| where (-A N )%w = v. 

k&z d \{ 0 } 

For (I74[) it thus suffices to show for arbitrary function w 

n | -t n ti 1 

||ry n A2U)|| < ||7y n+i VA2t(;||"+ T ||r(;||"+ T + ||tc||. 


By iterated application of Young’s inequality, it is easily seen that this family of interpolation es¬ 
timates indexed by even n follows from the following two-tier family of interpolation inequalities 
index by m € N 

||? ? 2m A m u;|| < ||r ? 2m+1 VA m u;||^||77 2m - 1 VA m - 1 u;||^ + || r? 2 ™- 1 V A m - 1 w|| 

and 

||r ? 2 m - 1 VA m - 1 HI < \\v 2 m A m w\\^\\r] 2 m - 2 A m - l w\\^ + \\r ] 2m - 2 A^w]]. 

Obviously, this two-tier family reduces to the two estimates 

||r? 2m Ar,|| < || r y 2m + 1 vA?;||^||r/ 2m_1 Vr;||^ + Hry 2 ^” 1 V^||, 

||f/ 2m_1 Vn|| < \\rj 2 m Av\\^\\r] 2 m ~ 2 v\\^ + ||r7 2m - 2 n||, 


which by Young’s inequality follow from 

h 2m M\ < (|| ? ? 2m+1 VAn|| + W^AvW^W^VvWA 
h 2m -^v\\ < {\W m Av\\ + \\r, 2m - l Wv\\) l 2 \\rj 2m -‘ 2 n||i 

Thanks to (ESI), these two last estimates immediately follow from integration by parts (the 
cut-off rj suppresses boundary terms), the Cauchy-Schwarz and the triangle inequalities. 
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